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New physics in the chromomagnetic flavor changing transition s — > dg* can avoid the strong GIM 
suppression of the Standard Model and lead to large contributions to CP-violating observables, 
in particular to the e' parameter, that we address here. We discuss the case of the Left-Right 
symmetric models, where this contribution implies bounds on the phases of the right-handed quark 
mixing matrix, or in generic models with large phases a strong bound on the Left-Right symmetry 
scale. To the leading order, a numeric formula for e' as a function of the short-distance coefficients 
for a wide class of models of new physics is given. 



Flavor-changing processes still offer one of the best means 
for spotting signs of physics beyond the Standard Model 
(SM) . The K decays into pions are among the best stud- 
ied channels both experimentally and theoretically, and 
despite the hadronic uncertainties in the theoretical pre- 
dictions, can serve as a great tool for probing new physics. 
The reason is that for a number of observables, in par- 
ticular the CP violating ones, the SM contribution is ex- 
tremely small, mainly due to the GIM mechanism [1 and 
the fact that the CP violating phase in the SM is sup- 
pressed by the smallness of the quark mixing angles. In 
turn, in the SM the GIM mechanism is intimately re- 
lated to the chiral nature of the weak interactions. As a 
result, probes of processes involving the GIM mechanism 
are well suited to test for nonchiral interactions. 

This is paradigmatic in one popular extension of the 
SM, Left-Right (LR) symmetry [2], which altogether 
gives a framework for restoration of parity in fundamen- 
tal interactions, nonzero neutrino masses, as well as vio- 
lation of lepton number and flavor [5] both at the reach of 
the coming round of experiments, fitting especially well 
with the scenario of TeV scale LR-symmetry [IHS]. The 
related direct searches at LHC, with important signa- 
tures through the new interactions and same-sign dilep- 
tons [7] , can explore this possibility up to ~6 TeV [5] and 
are already beginning to probe this interesting region 0|- 
[TT| . It is then important to assess the bounds on the 
model from existing phenomena. 

In the LR models, based the SU{2)l x SU{2)r x 
U{1)b-l gauge group, modifications of GIM are mainly 
due to the new right-handed gauge boson Wr and to 



its mixing with the standard weak gauge boson Wl- 
Bounds on the scale of the new right-handed gauge in- 
teraction were already addressed since the early days, a 
notorious example being the AMk box diagram [T^l |T3] 
where the GIM enhancement adds to a chiral enhance- 
ment of the matrix elements, and still leads today to the 
strongest bound on the scale of LR-symmetry, 
2.5-3 TeV 0]. Similar effects hold for the CP-violation 
parameter e jH]. The bottom line is that the inter- 
play of nonchiral interactions with the hierarchy of quark 
masses and mixings can lead to dramatic effects in loop 
diagrams. This is especially true in the phenomenol- 
ogy of strange mesons, and in particular for the direct 
CP-violation parameter e' and the chromomagnetic loop, 
that we address. 

The gluonic penguin operators have been traditionally 
associated to e'/e, because they pointed immediately to a 
possibly sizable effect. However, in the SM a partial can- 
cellation between the dominant gluonic and electroweak 
penguins translates into a large theoretical uncertainty, 
linked to hadronic matrix elements. For a review on 
the evaluation of e'/e and additional literature we refer 
to P^17j . In any case, e' is naturally tiny in the SM, 
and can serve as an efficient tool for constraining new 
physics. 

In this work, we address the contributions of nonchiral 
interactions in the chromomagnetic operator, and its ef- 
fect on e'. In the analysis we first give a parametrization 
of the effects of new physics in e' which is applicable to 
a wide class of models with nonchiral interactions. In 
the context of Left-Right symmetry, as is known current- 
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current operators mediated through the left-right gauge 
boson mixing gives large contributions to the K — >■ tttt 
amplitude. This issue was studied in detail in [il fTBH^ . 
together with the other flavor constraints on the model. 
However, the effect of the chromomagnetic operator was 
not considered. We study its impact due to the effec- 
tive K TTTT transition whose hadronic matrix element 
computed within the Chiral Quark Model (xQM) |21j. 

The xQ-M^ provides an interpolation between short- 
distance QCD and its effective description in terms of 
the octet of Goldstone mesons, below the scale of chi- 
ral symmetry breaking (for a recent discussion see |22)'). 
The chiral lagrangian coefficients are determined order 
by order in the momentum expansion by integration of 
the constituent quarks and depend on three non pertur- 
bative parameters: the constituent quark mass and the 
quark and gluon condensates. Via a fit of the A J =1/2 
rule va K ^ TTTT decays, the authors of ref. [15] obtained 
a non trivial phenomenological determination of these 
three parameters that allowed for a correlated calcula- 
tion of e'/e and of the AS' = 2 bag parameter Bk within 
the xQM approach, at next-to-leading order (NLO) in 
the chiral expansion p3Vl25j . We will use these values of 
the parameters in our analysis. 

For the LR model, we shall see that only the chro- 
momagnetic operator plays a dominant role in e', once 
other existing constraints from K and B physics are con- 
sidered, and implies a bound on the free phases involved, 
in the hypothesis of TeV scale LR-symmetry. For other 
new-physics models, constraints from e' via the chromo- 
magnetic operator were studied in [55H55] . 

The paper is organized as follows: In section |T] we 
describe the effective operators involved when nonchiral 
new physics is present, including the dipole ones. We also 
review the short distance coefficients in the case of the LR 
theory. In section [llj by running with the mixed anoma- 
lous dimensions, we compute the Wilson coefficients at 
0.8 GeV, which is our chosen scale for matching with chi- 



ral perturbation theory. In section III we describe (and 
update) the bosonization of the chromomagnetic opera- 
tor. This enables us to make contact with the K — >■ tttt 



amplitude and with e' in section IV in general and in 
section |V] for the LR model. In section IVII we draw our 
conclusions. 



the Left-Right model, the standard set of operators has 
to be extended. In the case of AS = 1 discussed here the 
complete set at low energy involves 28 operators. 
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The notation is (qq)L,R = qip,[^Tl5)q, L,R = (l=F75)/2, 
and the summation on q = u,d, s is implicit. f ^^'^ 
SM operators usually denoted as (5i,2- The dipole opera- 
tors Qg^-y are normalized with m^, for an easy comparison 
with existing calculations, and for keeping unaltered the 
anomalous dimension. It is known that some of the oper- 
ators above are accompanied by an enhancement due to 
the different chiral structure, either in the short distance 
coefficient, in the running, or in the matrix element. Such 
situation occurs for the gluo-magnetic operators Qg'^, as 
we describe below. 

At leading order the operators generated by the SM 
and the LR short distance physics are: Q^^, Qa, Q'4, 
Qe, Q'e, Q7, Q'7, Q9, Qj. with A,B = L,R. 

Their coefficients are: the current-current ones 
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I. NEW PHYSICS 



the penguin ones 



The effective lagrangian for flavor changing can be writ- 
ten in the form Las=i = —{GF/V^)J2i^iQi + ^-C-: 
where Qi are the relevant operators and Ci the corre- 
sponding coefficients (and Gp the Fermi constant). In 
the Standard Model the AS = 1 processes are usually 
described by a (over)complete set of operators [26l l27]. 
They involve tree-level operators as well as QED and 
QCD penguins. In models with both chiralities such as 
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and the dipole ones 



m,C^= E, [m^Af^Sa'-'' +Cm,Af^F2^^ +md/3Af ^Sa^^] 



(4) 



e„ = 2/3 is the u-quark charge, 



In the above 

mf/rn^^ with i ^ u,c, t, and -F^^-j and i?^^ are the loop 
functions, given in appendix A Then, f3 — My^^/Myy^ 
is the ratio of the electroweakto the LR scale and C, is 
the Wl-Wr mixing. Note that in the hypothesis of LR 
symmetry at TeV scale, /? ~ 10~^. Also C is of order 
/3 or less; for instance in the minimal LR models it is 
C ~ with X < 1 the (modulus of the) ra- 

tio of the two VEVs of the Higgs bi-doublet, and a its 
phase. We will consider below the specific case of min- 
imal LR models, referring to [4j [28] for definitions and 
details. Finally Xf^ = V*/V^^, where Vl and Vr are 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix and its 
right-handed analogue. A crucial new ingredient in Vr 
is the presence of (five) additional phases, besides the 
Dirac one. These can be parametrized as {U = u^c,t, 
D = d,s,b) 



''UD — 6 Vj 



R 
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with Vr the mixing matrix in standard CKM form. 

The terms in the expressions (H, (H, Q for the coef- 
ficients should be understood as generated at the decou- 
pling of the relevant heavy states, and thus at different 
scales, namely: or rrit for the current-current and 

top-dominated loops, rric for the charm dominated loops 
etc, and mwj^ for the RR current-current. 

A similar set of operators Qi2c with the c-quark re- 
placing u, is also generated by the short distance physics, 
and also by renormalization at scales larger than rric. 
On the other hand, the further operators involving the 
t quark are not explicitly required: for the LL and RR 
operators this is due to the GIM cancelation above nit 
(also in the running); for the LR ones, they are only gen- 
erated at electroweak scale through the LR-mixing 
Lastly, there are also penguin operators built through 
the LR-mixing, which are chirally suppressed and give 
subleading (negligible) contributions. 

^ Clearly, if Qf2t w^re generated at high scale, they should be 
taken into account, because due to the mixed chirality a GIM 
cancelation is not effective in the running (see also comments in 
appendix [b|| . Also, some more operators of the form (sq!)^ (Jd -|- 
Ss)l,r mixing with the penguins are present in general, see 29], 
but they are not generated in the LR model and it is also difficult 
to generate them in models where new physics sets in at scales 
higher than the electroweak scale. 



From Q it can be seen that the coefficients Cg'^ re- 
ceive a large contribution in the LR model, due to the 
different GIM mechanism. In fact, the mass insertion on 
the external fermion legs in the SM (rUs) is replaced in 
the LR model by a mass insertion inside the loop (to^). 
The loop is then dominated by nic leading to an enhance- 
ment of rric/ms ~ 100. In addition, the factor Xf^ /Xf^ 
gives a further large enhancement of 10"^ , which compen- 
sates the LR-scale suppression (. Both Qg and are 
present, and the LR contribution ends up being a factor 
^ 200 larger than the SM one, at short distance: 



2 X 10^ C - 200 . 



msFi^ixt)\V,*,Vt,\ 

(6) 

Here, the factor 2 accounts for the contributions LR+RL, 
and we considered Vr c:^Vl, which is a general prediction 
of minimal LR models 

The new phases ^ containted in Vr, together with 
the enhancement above, can directly induce a sizeable 
CP violation. It is therefore important to address the ef- 
fect of this operator on e', which we study along the lines 
of [301 EI] ■ In order to deal with this low energy phe- 
nomenon, two steps are necessary: the first is to renor- 
malize the coefficients at low energy, in the range of chiral 
perturbation theory; the second is to use the matrix el- 
ements (27r|(5i|ifo)j or equivalently to match with chiral 
perturbation theory. In the following section we renor- 
malize the full set of coefficients down to the scale of 
0.8 GeV, and in the next we match with the chiral la- 
grangian in the context of the Chiral Quark Model. 

The need to evaluate the Wilson coefficients at such a 
low QCD scale is dictated by the requirement to use the 
matrix elements calculated in the context of Chiral Quark 
Model in Chiral Perturbation Theory, whose cut-off is the 
chiral symmetry breaking scale. In order to assess quan- 
titatively the scale dependence of the result, we remark 
that by varying the matching scale between 0.8 and 1 
GeV the chromomagnetic Cg'^ and the current-current 
ones Ci2'^^ vary at most by 5% and 10%, respectively.. 
These uncertainties are well below those of the matrix 
elements. 

II. RUNNING TO LOW SCALE 

The mixing of operators ([T]) is described in detail in ap- 
pendix[Bj At leading order, the operators can be split into 
two sets, of opposite chiralities, corresponding to the two 
columns in ([T]). The low energy coefficients together with 
the matrix elements of all the operators are also sufficient 
to give an estimate of the impact on e' for quite a large 
class of models of new physics. This will be presented in 
section IIVI 

In the particular case of the LR model, the low en- 
ergy coefficients are shown in table [l] The running takes 
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TABLE I. Coefficients for the dominant new operators in the 
minimal LR model, evaluated at fj, — 0.8 GeV. 

into account the whole set of operators including the 
SM penguins, but we show the result only for the op- 
erators containing the LR scale /? or C which have an 
impact on e', and the magnetic operator Qi^'^ which 
is also enhanced. The results are normalized to A„, 
to compare with existing calculations. The coefhcients 
Cg'^ , compared with the complex part of the SM result, 
C^{SM) ~ 0.34Af PO], confirm the important role of Qg 
from new physics. 

In the detail of the running it is worth noting 
that, despite the reduction of ~0.5 due to their own 
anomalous dimension, C^'^ receive contributions from 
^LL,LRML,RR ^ rpj^^ largest additional contribution is due 

RL LR 

to ^1 2c scales above rric, while the contributions 

from C^2u^^ are suppressed by the u quark mass. This 
is due to the internal mass insertion in the (two) loop 
graphs responsible for the operator mixing in the anoma- 
lous dimension matrix, and is an other consequence of the 
nonchiral nature of these operators. As a side result, this 
additional contribution preserves the same combination 
of phases appearing in the original short-distance Cg'^ . 

From table [l] we can calculate the contributions to e' 
of these operators. This requires the evaluation of matrix 
elements which we review now for Qg '^. 

III. BOSONIZATION OF Qg 

The bosonization of Qg was addressed in [301 in the con- 
text of the Chiral Quark Model. Here we review the 
computation, which leads a minor numerical correction. 

Under chiral SU{3)l x SU{3)r rotations the Qg op- 
erators transform as (3^,3^), and thus they give rise 
to particular terms in the chiral Lagrangian. While by 
symmetry arguments there are diverse possibilities (see 
e.g. [32j in naive dimensional analysis) in the context of 
the Chiral Quark Model, only one form arises [50] • This 
is true in the SM as with the separate operators Qg, Q^. 
One has 

Cq^ = Tr [(EtxA_ + A_AtE) D^J^^D^Y.] , (7) 

where A_ = (Ag — iXj)/2, and where the matrix of 
two coefhcients X = diag(0, G|^, Gg ) replaces the sin- 




FIG. 1. Diagrams for the bosonization of Qg, in the fixed 
point gauge. Note the flow of gluon momentum. 

gle coefficient G^"^ and the running quark mass matrix 
M — diag(TO„, md, TOs) of the analogous calculation in 
the SM. In fact, the coefficients C^^--^ of the AS = I 
transition induce a tiny breaking of the chiral symmetry 
which plays the same role as the s and d quark masses in 
the SM. We also observe that from the point of view of 
chiral perturbation theory the above lagrangian is 0{p'^) 
in the LL and RR terms proportional to light quark 
masses, but 0{p^) in the LR and RL terms, proportional 
to mc- In any case, the coefficients Gg'^ will respec- 
tively be proportional to Cg'^. To determine them, it is 
convenient to evaluate some amplitude both through the 
above chiral lagrangian and in the Chiral Quark Model, 
and compare the results [30]. We do this in the 'unro- 
tated' picture [T7]. The simplest process is the off-shell 
— > TT^ transition, which at one loop is given by the 
three diagrams shown in figure [T] The two external gluon 
lines are attached to the giuons in the external thermal 
QCD vacuum and lead to a coefficient proportional to the 
gluon condensate. To deal with the thermal and color 
average, it is best to adopt Fock-Schwinger fixed-point 
gauge {x^'Afj, = 0) [53]. Due to this gauge choice transla- 
tions are broken and two fixed 'sink' points for the gluon 
momentum are defined, chosen here to be a; = and x 
at the K and tt vertices. Then the three diagrams in 
figure [T] are 

a =tr[(^ - + m)4i^ + m)qg{^ + m)(^ + to)] 
A[(7, m]^A[(7 — k, to] A[g — p, to] 

b =ti-[(i - ^ + m)qg{i + m)4{4 + to)(^ - / + to)] 
A[(7, to] A[g — k, m]'^A[q — p, to] 

c =tr[(^ + m)qg{^ + to)(^ + m)4{4 + ^ - ^ + m)] 
A[(7, to]^A[(7 — p, to] A[(7 — p -I- fc, to], (8) 

where p is the K and tt momentum; k is the in- 
coming gluon momentum; A[p, to] = ~ to^); 
qg = -iGF9snis{C^L + R)G°'^ a ^fi 1 4^ and to is 
the constituent quark mass. Also (see [SB]) 4 = 
{gs/2)jf^Gfi^d/dki, where a derivative with respect to 
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the gluon momentum has to be taken, after which k is 
set to zero|^ By the same prescriptions one shows that 
no external gluon momentum flows through the chromo- 
magnetic operator. Finally, the two gluon field-strengths 
are averaged in the gluon condensate, Ga^pG^^s 
(7rV6gs) (^GG) {ga^gps - gasgfs-i)- 

The loop integration produces a term at order zero 
in p which is canceled in the leading chiral lagrangian 
in agreement with the FKW theorem [351. The second 
order term gives the desired result 



iGf 



a 



,{7 +1+8) 



(9) 

the three numbers being relative to diagrams a, 5, c. 
Comparing this amplitude with the one calculated from 
the chiral lagrangian ([t]) , one finally finds the coefficients 



12m \ TT 



y/2 167r2 



(10) 



The factor of 2 corrects 11/4 appearing in [3U], and leads 
to a 30% reduction of the matrix element. This result 
is also confirmed by a similar calculation within the 'ro- 
tated' picture. 

Before using this result for the calculation of e', let 
us note that the additional contribution from the off- 
shell chromomagnetic sdg vertex, shown in [3^ to be of 
the same order as Qg in the SM, is strongly suppressed 
in the case of nonchiral interactions (again because the 
mass insertion happens inside the W loop). 

Finally, a double insertion of Qg, leads directly to 
(long-distance) AS* = 2 processes and can be calculated 
similarly. However, the process is doubly loop and Gp 
suppressed and the result is negligible for both AMk and 
e (see appendix [Pj). 

IV. RESULT FOR e' 

The direct CP parameter e' is defined by 



i / ImA2 



ImAo 

RbAq j p 



(11) 



where p, q are the -fC", K mixing parameters and 
Lo = A2/A0 ~ 1/22.2. The ratio p/q ~ 1 with an excel- 
lent approximation. The isospin amplitudes Aj (/ = 0, 2) 
are defined from the AS* = 1 effective Hamiltonian as 
((27r)/|(-i)i7AS=i|i^'') = Aje'^', where Si are the strong 
phases of tttt scattering. We calculate the imaginary 



^ Due to the absence of translational invariance, the use of the 
one-external-gluon effective quark propagator in the fixed point 
gauge I37| is not correct in diagram (b), where the gluon mo- 
mentum flows to the origin {K) passing through the operator 
insertion. This leads to a mismatch between diagrams a and b 
and to a numerical correction of the result in ref. I30| . 



part of the amplitudes, while for the real part we take 
the experimental value: ReAo = 3.33 x 10~^ GeV and 
ReA2 = 1.49 X 10"^ GeV. 

The amplitudes Aq and A2 for the standard operators 
are collected in appendix Icl The ones of Qg'^ are easily 
calculated from the chiraTlagrangian ([t]). One has the 
isospin decomposition 



a: 



-At^'"' 



a: 



1 0^+^ 

, (12) 



where the amplitude for K'^ — > tt+tt is 
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6mf^ 
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(13) 



and the one for K'^ — > tt^tt'^ vanishes. In the follow- 
ing we use / = 93MeV, and for the gluon condensate 
and constituent quark mass we adopt the central values 
{^GG) = (334MeV)^ and m = 200 MeV [H], obtained 
by consistently fitting in the model the A J =1/2 selec- 
tion rule. 

Using the running and the matrix elements, one can 
generically describe the contributions to e' of the differ- 
ent operators, as a weighted sum of the coefficients con- 
tributing from the desired scales 



e = 



^ ^ Wi (Mn) Im Gi(/z„) 



(14) 



where, at each scale /j,„, the Gi(/i„) are the coefficients 
and Wilpn) their weights. The scales fj,„ can be either 
taken as the ones where the short distance coefficients 
are generated, i.e. mw, rnc, etc., in which case the Wi 
account for the running and the matrix elements, or some 
low energy scale if one includes the running in the G^. In 
table [n] we collect the numeric values of Wi computed by 
taking into account the complete running from a choice 
of different scales down to 0.8 GeV, together with the 
required matrix elements (see appendix [C|) . 

As discussed in appendixjC] the hadronic uncertainties 
present in the B factors of departure from vacuum satu- 
ration approximation can be sizable and may vary from 
10 to 50% for the better known operators, to a factor of 
order one for Qg and Qff . For the SM operators the val- 
ues adopted in table |ll] are taken from the Chiral Quark 
Model calculation [15 . For the LR operators Q^f", a 
determination is missing but an guess can be given by 
noting their similarity with operators Qt^s- 

For Qg'^, while the leading order chiral bosoniza- 



tion ([7]) results in a m^/m]^ suppression (see (13)), this 
may cease to be true in higher orders and may result in a 
further enhancement. Together with chiral loops, this is 
likely to lead to order one correction coefficients Bq 2 to 
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0. 


0. 
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15. 


0. 


0. 


1.9 


1.0 


0.16 


-0.74 


2.7 


8.0 


32. 


110. 


2.8 


1.9 


0.038 




= 0.8 GeV 


0. 


0. 


42. 


11. 


0. 


0. 


2.0 


1.2 


0.22 


-0.55 


2.2 


6.7 


25. 


88. 


2.8 


2.1 


0.043 



TABLE II. Weigths Wi (times 10) of the coefficients d in the determination of e', applicable to coefficients from different scales. 
The weights relative to the opposite chirality operators have opposite sign. 



be added to eqs. ( 12 ). While we stress the need for a ded- 



icated assessment of these corrections in the view of new 
physics, in the present analysis we conservatively assume 
— 1 J keeping in mind that a possible enhancement 
would make our bounds below stronger. 

Clearly, taking the central values of the hadronic ma- 
trix elements is sufficient for the scope of assessing the 
relevance of contributions beyond the SM, and, in partic- 
ular, in view of the leading role of the chromomagnetic 
operators shown by the present analysis. Also, the renor- 
malization evolution is performed at leading order, as is 
the determination of the starting conditions in eqs. ([2]- 
[4]). Nevertheless, since the penguins can be considered 
as NLO contributions, a NLO correction to the current- 
current starting conditions was also inserted (see j27|). 
In this respect let us again remark that while a NLO 
analysis is necessary for the SM, it is not crucial for as- 
sessing the constraint from e' on new physics. We believe 
table [n] with formula ( 14 ) to be useful in analyzing the 
impact of e' for quite a wide class of new physics models. 

V. CONSEQUENCES FOR THE LR MODEL 

We can finally estimate the numerical impact on e' of 
the new physics operators in the case of LR-symmetry, 
using the values of the low energy coefficients summarized 
in table |l] and the last line of table [ll] (or equivalently 
the appropriate short distance coefficients with the first 
lines). We find 



|C|1.25[sin(a-e„-0<j)-f sin(a-0„-0,)] 
+ Id 0.010[sin(a -9^ -9^) + sin(a -0^ ~0,)] 



/3 0.013 sin(6ld -Os) 



(15) 



where the first line is due to the dominant Qi^'^^ , the 
second to Q^'''^, and the last to Qf^f The penguin 
contributions are only responsible for minor corrections 
in the above numeric result. 

For TeV LR-scale (/3 ~ 10^'^), we see that the above 
contributions can give overdominant contributions to e'. 



even having assumed similar left and right quark mixing 
angles, as in the minimal LR models. This is true both 
for Qg'^ and for the other operators. Their impact may 
be different and it depends, in addition to the LR scale 
/3, on the CP phase a and on the extra phases in Vr. The 
actual implications for a given LR model depend thus on 
the available freedom in choosing these phases. Let us 
describe the bounds in different scenarios, assuming that 
the new physics can contribute as much as 100% to e'. 

Assuming generic 0(1) free phases, the LR scale is 
constrained to lie above a large limit Mw^ ^ 25 TeV. In 
fact this amounts to the highest limit on the right-handed 
scale. It is also worth recalling that for order one phases 
another large bound of about Mwr ^ 15TeV results from 
e, while limits from B mass difference and CP violation 
are less stringent P|, . The argument can however be 



turned around and ( 15 1 can be used to put constraints 



on the phases, in the scenario of TeV LR-symmetry. 

In the minimal LR models, phases are either strictly 
predicted or are free, depending on the choice of LR- 
symmetry, which can be generalized parity (P, exchang- 
ing fermions ipL ^ V'i?) or generalized charge conjugation 
(C, exchanging fermions V'l V-'flJi following the analy- 
sis of ref. 

An important common constraint resulting in both 
cases from e is that 9d — Og is close to (or marginally 
to tt) at least as 10^^. This implies that the contribu- 
tion of Qf f' to e' can be neglected (the last line in (15)). 



This can already be seen from the low energy coefficients 
in table |l) From the same table |l] it is also important to 
note an other consequence of the e constraint, namely 
that thanks to dg — 9d the coefficients are practically 
complex conjugated under L ^ R exchange. As a re- 
sult, the contributions to e', which are proportional to 
L 



"9 



C^, see 



13 ) are purely 
On 



^ Also A^^ is complex, but its phase is 0(1/1000) and is subleading 
in this expression. 



R combinations (as C, 
imaginary and thus with maximal imaginary part 
the contrary, only the combinations L + R would enter in 
other CP-violating observables involving an even num- 
ber of mesons, making the imaginary part suppressed by 
&s — Od S lO^'^. This is the case for instance for the 
contribution of the magnetic operators to the CP 

asymmetry in if — > 7re"'"e~, which is thus suppressed, 
despite the enhancement of the Wilson coefficient. This 
situation can be contrasted with the one occurring e.g. in 
supersymmetric models, where a correlation between e' 
and K — ^ 7r£+£~ can be inferred ESI- Nevertheless, the 
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enhancement of C^'^ would survive in CP-asymmetries 
with an odd number of mesons, hke K 7nr£^£~ (whose 
analysis brings in the leptonic sector of the LR models 
and is beyond our study). 

Let us then describe the impact of the first two lines 



the second line, due to Qg For e'^^ to at most saturate 
the experimental value one has the constraint 



in ( 15 1 to e', for the two possible choices of LR-symmetry. 



In the case of V, due to the hermiticity of the Yukawa 
couplings, the phases in Vr are all predicted in terms 
of the phase a, and they are all close to or tt. The 
neutron EDM then poses a strong constraint which to- 
gether with e and e' leads to the strong limit Mw^ > 
8-10 TeV [in]. As discussed in g], a TeV-scale LR sym- 
metry is still allowed by resorting to an unappealing fine- 
tuning with the QCD strong phase 9. In this case, the 
e' gives alone a bound, because from e one must have 
X sm a ~ 

10-3 gj Hn], so that the hmit x i.e. C ^ 



that would suppress e' in ( 15 1 is not permitted. Then, by 



using the values for the predicted phases (see eq. (29) in 
rcf. [19 ) and exploiting conservatively the free signs, in 
particular u and c opposite, one finds the numeric resull[^ 

|e'^^|^5.710-6(/3/io-3). (16) 

This result holds in the natural regime x < 0.1, where 
analytic expressions for the phases are available. Com- 
paring e^^ with 100% (50%) of \e'\exp. - 3.92 x lO^^, we 
obtain the constraint Mw^ ^ 3 (4) TeV. Here the main 
contribution comes from Qi^'^^ , while Qg'^ contribute 
subdominantly, with the effect of softening the limit with 
respect to the conclusions of ^4jj^ 

In the more interesting case of C as LR-symmetry, the 
phases ^ in are free. This time, the bound from e to- 
gether with the Bd,s systems put the stronger constraint 
9d — Os < 10^^ 4 . As a result, from (15) one has 



Id 2.5Osm{a-0u-0d) 



+0.020sin{a-ec-0d) . (17) 

The first line due to Q^2'^^ is dominant, but one can 
note that the phase combination appearing there is inde- 
pendently constrained by the neutron EDM: as studied 
in ;4l Hg ED] one has |C|sin(a-6l„ -Od) < 10"^. This 
implies that the current-current contribution can be ne- 
glected here, and thus we are left with the dominance of 



This corresponds to the case — 6s| ~ 10~^. There is also a sec- 
ond possibiUty, with — 9s| ~ 7r + O(10— and xsina 10"^, 
but it is disfavored by the B^j „ mass differences and CP-violation. 
Moreover, in this case there are cancelations in each line of ( |15| l 
and the situation is more ambiguous, since it depends on the 
precision to which the equality ~ Vn of mixing angles holds. 
We recall that with a global numerical study of this model T, 
the angles could be deviated as much as 20%, which would spoil 
these cancelations in e' . Therefore also in this case one can ex- 
pect a dominant contribution as in llml. 
^ We also do not agree with the strong bound derived in ref. |20| . 
using only the isospin-2 amplitude of the operator Qf^. 



\C\\sm{a-e,-ed)\ <2.0x 10" 



(18) 



(or a correspondingly more stringent one for a subdomi- 
nant |e'|L/j). This represents a correlated bound between 
the phases and the LR-symmetry scale/mixing. For un- 
constrained phases 6c + Od one would require Mw,, 
8 TeV, or vanishing LR-mixing, while for TeV-scale LR- 
symmetry and large LR mixing one puts a constraint of 
~ 10^^ on the involved combination of phases. 

This constraint is relatively stronger than the one re- 
ported in ref. jH] , where a different evaluation of the Qg 
matrix element was adopted, and where only the 'charm' 
couplings were considered. As discussed above, the un- 
certainty (possible enhancement) in the matrix element 
of Qg could make this bound even stronger by an order 
one factor. 



The bound ( 18 ) is also stronger than the analogous 
limit inferred from the s — >■ c?7 decays [411 02], which 
through Qi^'^ involve the same enhancement and the 
same phases as the chromoniagnetic operator. 

VI. SUMMARY AND CONCLUSIONS 

In this work we addressed the effect on e' of new physics 
in the chromomagnetic dipole operators, which can have 
a huge enhancement with respect to the SM, especially in 
the presence of nonchiral interactions. The paradigmatic 
example for this effect appears in minimal Left-Right 
symmetric theories, where the Wl-Wr gauge-boson mix- 
ing leads to an enhancement of 10^ in the short distance 
loop coefhcient, so that even with a scale of new physics 
in the TeV region, an enhancement of two orders of mag- 
nitude results. Together with the presence of new phases 
in the Right quark mixing matrix, this can lead to a dra- 
matic impact in e'. 

To evaluate quantitatively the effect, we considered the 
dipole operators together with the full set of four quark 
operators which can give rise to CP violation in _ftr — > tttt 
decays. We considered their renormalization and mixing 
(at leading order) from short distance to the low scale of 
matching with chiral perturbation theory, where the ma- 
trix elements can be estimated. For the chromomagnetic 
dipole operators we reevaluated the corresponding ma- 
trix element in the context of the Chiral Quark Model 
(correcting the previous existing calculation). Tor the 
SM operators we adopted the estimates consistently de- 
termined in previous analysis of e'/e |15) . These were also 
used for an estimate for the LR current-current operators 
Qf f (see discussion in section [iv| and appendix [cj). The 
set of high energy operators considered is fairly complete, 
and can serve also for estimates of the impact of other 
new physics models on e'. 
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We applied the results to the case of the minimal Left- 
Right model, showing that e' receives contributions from 
the chromomagnetic operator as well as from the current- 
current ones Q\^'^^ . These are in general large, but we 
noted that they are severely constrained by the nEDM, 
with the result that the chromomagnetic operators turn 
out to be dominant. One can expect this to be a fairly 
generic situation, because new CP phases contributing to 
Qi 2 ' are usually contributing as well to the nEDM. 

In the LR model, focusing on the case of general- 
ized charge-conjugation C taken as LR-symmetry, where 
new phases are free, this allows us to derive a corre- 
lated constraint between the LR-gauge boson mixing (or 
LR-symmetry scale) and the relevant phases. The bound 
amounts to C ^ 10"'' (or Mwpi > 8 TeV) for arbitrary 
phases or, in the scenario of TeV-scale LR-symmetry, a 
constraint of 10~^'^ on the involved phases. 

The dominance of the Qg contribution can then lead 
to constraints on processes involving the same combina- 
tion of phases. This is true for instance for the mag- 
netic operators Oil'^ which are also GIM-enhanccd, and 



as we argued will enter in K ^ 7r7re+e~ CP- violating 
asymmetries. Also, (a — 9c — dd) enters in the charmed 
mesons physics the analysis of which is beyond the scope 
of this work. Nevertheless let us point out that it enters 
the decays of the D meson via c — )■ wj, whose short 
distance contribution is overwhelmed by the long dis- 
tance ones [43], but also it enters the CP-violation in 
the D — >■ KK, TTTT channels, for which anomalous signals 
have been reported by the LHCb collaboration [44 . The 
interesting analysis of the related charm physics in the 
LR-models will be the subject of a separate work. 
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Appendix A: Loop functions 



The loop functions relevant for the SM and the LR model are [TU . 

a;, (-18 + llx, (4 - 16a;, + 9x2) In x^^Sx? - 2a;,- - 6) 19a;3 - 25x2 4 

■F^i = — 7 -r. — — - 7^ rr , = — ^ Ina;, + — / + - Inx, 

^ 12(xj - 1)3 6(xj - 1)4 ' 1 18(xj - 1)4 ^ m{x.j - 1)3 9 ^ 

4(x,-l)4 ' 2 12(x,-l)3 2(x,-l)4 ''''^ 

-4 + 3x,-+x3-6x, lnx,- Sx^ - 31x, + 20 x,-(2-3x,) 

^2 = 2(^r3Tp ' - 6(x,-i)2 

where Xj = (y^)^, j ~ u,c,t. In addition to these, one has F(^2 = ^12(^^1) and similarly for 



Appendix B: Running of all AS — 1 operators 



For our purposes, the relevant operators are the Qi appearing in M plus the eight Qi2c where u quark is replaced 



by c. At leading order (LO) the LR operators mix only with Qg^ (m addition to themselves). Similarly, the RL ones 
mix only with Q^^- The operators can thus be divided in two decoupled sets of 18 operators each, related by the 
exchange L ^ R: 

Q^,Q!;}, (Bl) 

{QlcjQ2ctQl tQ2 tQictQ2c tQi : Q2 1 Qs' Q 4' Q 5' Q 6' Qt ' Q St Qd, Q 10' 

Q'\Q!^}, (B2) 



as in the two columns of eq. (|lj). The corresponding vectors of coefficients in the two sets, Cl,r{i^), evolve separately 
according to the renormalization group equation 



9 In /i dg ' 9 In 
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with i = u,s, c. The 18 x 18 anomalous dimension matrix 7 is the same in the L and R sectors and reads [171 US] 
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(B4) 



where n f and ri„ are respectively the number of active quarks, and of active up-type quarks. The mixing of the Qi . .^10 
operators with the (chromo) magnetic operators appears at NLO j48| . and we adopt the anomalous dimensions in the 
HV scheme. Note, the off-diagonal terms in the last two columns carry explicitly the ratio of mass insertions responsible 
of the operator mixing. In fact, while for the mixing of the LL and penguin operators with the magnetic ones the 
mass insertion on the external legs is and coincides with the normalization of Qg,-y, for the mixing Q^^'^^ — )> Qg.^ 
the mass insertion is that of the internal quark (m„ or TOc), breaking the usual GIM cancelations. It also follows that 
the traditional description in terms of the and Zi variables |27j is no longer appropriate, and we need to perform 
the running of the whole vector of coefficients Cl (and Cr). Clearly, for the SM part the results coincide, as the GIM 
cancelation is effective in the mixing with the dipole operators. Finally, the operators involving the charm quark are 
integrated out at their threshold, and accordingly the anomalous dimension matrix is projected, below this scale, on 
the remaining set of low energy operators Q^f , Q3-10, Qg, Qj- 

We perform the running choosing as{Mz) — 0.1176, and with starting coefficients introduced separately at the 
relative scales of decoupling. The running is performed down to 0.8 GeV where the matrix elements are evaluated. 
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Appendix C: Amplitudes 

For all the operators the amplitudes Aq^2 are expressed in terms of their K — )■ (7r7r)/=o,2 matrix elements {Qi)o.2' 

^=^Q(Q,)o, A2=Y,C^{Q^)2■ (CI) 

i i 

We report here the matrix elements for the relevant q^^-^^^^^^^^ [15] : 

3v6 3v3 

(Q")o = -(QD2= ^^^0, (Q")2=-(0n2 = ^^S2' (C2) 



)f ^>o = -(Qf - ^^'^3^ + '^^ , (Qf ^>2 = -(Qr>2 = ^ VS - 6^) ■ 



with 



(Q2'^>o = -(QDo = ^ 3^ ^ So^-"" , (Qr')2 = -(Qr>2 = 3 V 3 - 2Z) , (C3) 



X = -(7r-|d7^75u|0)(7r+|?27''s|i4'°) = i^/2f^{m\ - ml) ~ O.OSi GeV^ 
Y = -(7r+7r^|uu|0)(0|J75s|K°) = iV^fKA^ ~ 0.22iGeV^ 



Z = -(7r~|d75u|0)(7r+|?Is|A' ) = i\/2/^A^ ~ 0.18i GeV^ (C4) 

where A = m\- / [ms + md), and /7r,/f the tt and K decay constants, and the quark masses are evaluated at ^ = 0.8 GeV 
(i.e. rus ~ 200 MeV). Since we use the matrix elements at /i = 0.8 GeV, also the Bi coefficients of departure from 
vacuum saturation have to be evaluated at this scale. 

The SM ones, determined in the Chiral Quark Model via a phenomenological approach based on the fit of the 
A/ =1/2 rule in iC — tttt decays, can be taken from ref. [T5] (see table VI), where one can also find the "correlated" 
matrix elements for the operators (33,...,io- For the above current-current operators one finds the central values 
~ 9.5, Sq ~ 2.9, B2 — 0.41. For the gluonic and electromagnetic penguins relevant to e'/e it is found ~ 1.6 
and Bl ~ 0.92. 

Concerning the B^'^'^^, their evaluation is still lacking both in the Chiral Quark Model as on the lattice, and to 
our knowledge also in 1/Nc expansion. Some hints can be derived from the observation that the electromagnetic 
penguins Qj^s transform as (8l,8/j) as do the Qi2- Then, their leading bosonization and chiral loops coincide 
(see [23]) so that one can expect the B parameters of Qff to be very similar to those of Qsj- For the isospin-2 
amplitudes this correspondence has even been argued to be exact [20] so that using the results reported in [T5] we 

1 2 LR 7 S T 8 12 LR 

can set B2 ' — B2 — 0.92. For isospin-0 amplitudes, the larger Bq' ~ 2.5 hint for Bq' ' also larger than one. In 
general, this is in accordance with the strong phases from pion rescattering in final state interactions which point to a 
correction factor of ~ 1.4 [IHIIS^, and also more simply with the correction factor traditionally applied to Y in vacuum 
saturation to account for the renormalization to the K scale in the pion matrix element, (1 + rnj^ / A'^prp) ^ 1.5, which 



enhances the isospin-0 amplitudes in eq. ( C3 1 . Therefore for the present analysis we adopt the conservative choice of 
central values Bq'^'^^ ~ 2 with 0(1) uncertainty. For the Left-Right model, the impact of the isospin-0 amplitudes 



of Q12 is fortunately limited. The result in the first line of (15) changes by 20% within a 1-3 range of i?Q'^'^^. 



11 



Appendix D: Double insertion of Qg leading to AS — 2 



The gluonic dipole operators can give effects also on the AS = 2 processes, via its double insertion or via its insertion 
together with other AS = 1 operators, most notably the SM current-current operators Q12 which have large (real) 
coefficients. These combinations constitute true long-distance contributions. A first estimate was given in |51) and 
the resulting constraint is not relevant for AMx and is marginal for e [IT]. In the presence of the whole set of new 
physics operators also chiral loops with multiple possible insertions should be 
evaluated, and we leave that for a future analysis. However, a simple double 
insertion of Qg can be readily estimated. It leads to K-K mixing through the 
diagram on the right. The two external gluons are averaged in the vacuum 
ghion condensate, and the total mixing hamiltionian is: 



H 



KK 



8 Af|- G\ 2 
-m: 



3 p 



167r2 



K 



GG)KK 



(Dl) 



Considering from tablejljthe low energy values C*^'^ — 2.7|C|e**'"~^'^ ' we 
find the impact on the KK mixing to be negligibly small, despite the huge 
enhancement of the dipole loop: from AMk — Re{Hj^j^)/2MK we have 




AM 



< 10-^^ GeV , (for C < 10 



-3^ 



(D2) 



which is six orders of magnitude less than the experimental value. Similarly, also the effect on e is negligible: we have 
e ^ 0.3|Cp cos{2a){6d — 9s)- Since ( < 10^^ and 9d — Og is at most 10~^, this gives no constraint. 
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